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Abstract 

A gauge theory of second order in the derivatives of the auxiliary field is constructed 
following Utiyama's program. A novel field strength G = dF + fAF arises besides 
the one of the first order treatment, F = dA — OA + fAA. The associated conserved 
current is obtained. It has a new feature: topological terms are determined from 
local invariance requirements. Podolsky Generalized Eletrodynamics is derived as 
a particular case in which the Lagrangian of the gauge field is Lp oc G 2 . In this 
application the photon mass is estimated. The SU (N) infrared regime is analysed 
by means of Alekseev-Arbuzov-Baikov's Lagrangian. 
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1 Introduction 



Symmetry is one of the most important concepts in physics. Group theoretical 
methods are used to find complete solutions for many systems, particles were 
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found by symmetry arguments, from Noether's work we learn how to relate 
symmetries to conservation laws. In this context a natural question is: by 
imposing to a system a generalized version of some previous existent symmetry, 
would it be possible to determine new features? In particular, starting from 
a kinetic Lagrangian with a global symmetry, would it be possible generate 
interactions via localization of the symmetry? In other words, could local 
kinematical symmetries also imply dynamics? 

The general answer to these questions is yes, if the local symmetry is rep- 
resented by a general Lie group, and it was found by Utiyama [1] in 1956 
showing that it is possible to keep the action invariant under point dependent 
transformations of the matter field if one introduces a new field, the gauge 
potential, and derives a minimal coupling prescription which determines the 
fundamental interactions. Some years later, Ogievetski and Polubarinov [2] 
presented a criticism to this gauge principle by means of what is known to- 
day as B-Field formalism [3|H] . Notwithstanding, gauge principle remains as 
a cornerstone of modern physics. 

On the other hand, if one assumes that the equations of motion are derivable 
from some Lagrangian, a natural way to generalize a given theory is to suppose 
that the Lagrangian of the field contains terms involving derivatives of higher 
order. Poincare in 1901, when discussing the law of inertia, already called our 
attention to the importance of higher order equations [5J: "(■■■) The law of 
inertia (...) is not imposed on us a priori (■■■)■ If a body is not acted upon 
by a force, instead of supposing that its velocity is unchanged we may suppose 
that its position or its acceleration is unchanged. (...) in the second case, [we 
may suppose] that the variation of the acceleration of a body depends only on 
the position of the body and of neighbouring bodies, on their velocities and 
accelerations; or, in mathematical terms, the differential equations of motion 
would be of (...) third order (■■■)■" The non-singular Hamilton-Lagrange the- 
ory was extended to arbitrary order by Ostrogradsky [B] in 1850, generalizing 
the form of canonical momenta. Following this reasoning, Bopp [7] and Podol- 
sky [8] independently proposed a generalization of electrodynamics containing 
second order derivatives. Quantization of the theory resulted in finite energies 
in 1-loop approximation. This leads to the idea that some quantum field cor- 
rections would be simulated by the Podolsky's effective term. This generalized 
electrodynamics was able also to explain the 4/3 factor in Abrahaam-Lorentz 
theory [9] as an electromagnetic mass term by means of an appropriate gauge 
choice and a Lorentz invariant cutoff induced by quantum effects. This cutoff 
provides an natural mechanism to estimate the photon mass in the Podolsky 
generalized electrodynamics, as we will show in the section I5T21 - the presence 
of this massive mode does not violate the gauge invariance as it occurs in the 
first order approach, and it is an intrinsic feature of the second order gauge 
theory. Inspired by Podolsky's work, Green [TDj included another term involv- 
ing the second derivative obtaining a generalized meson-field theory with finite 
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energies at 1-loop. The relative success of these achievements motivated some 
authors to propose finite extensions of Quantum Chromodynamics (QCD) 
[TT] and also to advocate that higher order terms would be able to explain the 
quark confinement. The undesired feature of the Podolsky's theory is the lack 
of unitarity at 1-loop for its quantum version, and it is an open question if 
this characteristic can be ruled out in a nonperturbative scheme. Besides all 
these motivations we emphasize that, from a theoretical point of view, higher 
order theories have many interesting features that justify their study by itself. 

The scope of the present work is to construct a gauge theory for Lagrangians of 
second order on the auxiliary field, following Utiyama's procedure. Section [2] is 
devoted to fixing the notation, to reviewing how the gauge potential A appears 
and why it ensures the invariance of the action integral. In the next section, we 
assume a Lagrangian of the type L (A, OA, d 2 A) for the gauge potential and 
show that this Lagrangian must depend on A and its derivatives only through 
the usual field strength F = dA — OA + f AA and a second field strength 
G = OF + fAF, a new quantity that naturally turns up in the second order 
theory. The covariance of F and G under the gauge transformation is proved 
in the same section [3] and the Bianchi identities for these two objects are also 
deduced. 

After settling the basis of the general theory, we proceed to the analysis of 
the current derived from the total Lagrangian - the matter one plus the gauge 
Lagrangian - in section HI If Utiyama's definition of current is kept, one obtains 
a quasi conserved current instead of a conserved one. An alternative choice is 
to define the current in a different fashion to enforce conservation. 

As an application, in section I5.1[ it is demonstrated that the Podolsky Gen- 
eralized Electrodynamics can be derived from the second order gauge theory 
from a Lagrangian of the type L P oc G 2 . Section deals with a particular 
non-abelian case. 



2 Local Invariance and the Gauge Field 

Let Q A (x), (A = 1, 2, N) be a general matter field whose Lagrangian density 
is 



HQ 




dQ A 



dx^ ' 



with equation of motion given by 



dL 



-d, 



dL 



(1) 



dQ A 



0. 
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We postulate the action integral, 

L d 4 x , 

to be invariant under the global infinitesimal transformation: 

Q A ^Q A + 5Q A , 
SQ A = T A B e a Q B , (2) 

where e a is an infinitesimal parameter independent of x (a = 1,2, ...,n) and 
T a A B are constant matrices. In what follows, we assume that the transformation 
(T5]) belongs to a Lie group G dependent on n parameters e a . Then the structure 
constants f^ c are defined by 

[T a , Tb) A B = fat,T c A B) (3) 

satisfying 

f m f 1 _i_ f m f 1 _i_ f m f 1 = n (A) 

JabJmc'JbcJma'JcaJmb w ' V / 

which is the same as Jacobi identity, and 

fab = ~fb a' (5) 

in accordance with ([3]). 

From the invariance of the action under (T5]) in any spacetime volume Q, it 
results 

Using the independence of the parameters we find, 

In order to write (J7J) we consider a variation 5 which is strictly functional, i.e., 
the spacetime point is not changed. 

Rewriting ([6]) by using Leibniz rule, 

The field equation (CQ) must be satisfied, so that 

BL 
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Now let us consider the following infinitesimal transformation with a point 
dependent parameter e a (x) (a = 1,2, ...,n): 

5Q A (x) = T a A B e a (x)Q B , 

T A B = constant coefficients, (9) 



In this case, 



^^'W' < 10 > 



or 



We see in this case that 8L does not vanish - still assuming that (171) is valid 
even when e a = e a (x). 

If one wants to preserve invariance of the Lagrangian under (Q, it is necessary 
to introduce a new field pQ , called gauge potential, A a ^ (x) which transforms 

as 

M; = / a c //(x) + 9/, (12) 
appearing in a new Lagrangian L'(Q A , d^Q A , A a ) through the combination 

V,Q A = d,Q A - T a A B Q B A a . (13) 



This new object is covariant, since it transforms exactly as the original field, 
i.e., 

5{y,Q A )=T A B V,Q B e a {x), 
and it substitutes the ordinary derivative in the original Lagrangian, in a 
prescription named minimal coupling for description of the interaction, 

L'(Q A ,d tl Q A ,A a fl )=L(Q A ,V,Q A ). 

Usual derivatives are replaced, in the original Lagrangian, by the covariant 
derivative given by (|T3l) . 



3 Gauge Field Second Order Lagrangian 



Let us take a Lagrangian for the auxiliary field as a functional kernel dependent 
up to second order derivatives, 



LoiA'dyA'daduA*), d v A a 



dA a 



i-' 



dx L 
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We will impose that this kernel is invariant under (fT2|) . so 

9L 



°-^5 (d a d v A\) = 0. 



Substituting the transformation law for the gauge potential in this equation 
and since e a and its derivatives must be functionally independent, we are led 
to 

9L ° /X + + TT^TtWX - 0, (14) 



dA a " " ' " 



d(d v A^) chv » d(d a d v A^ 



(15) 



+ ^ ^^L/ aa + Q / Q o z: \ fcb A n = Q) (is) 



<9Ln <9Ln dLt) , , 

H 1 = 0. (17) 

' a /a a An \ ' a / a a An \ \ > 



d (d a d u A^) d {drdvA'a) d {d„d a A% 



This set of equations forms a hierarchy informing us about the dependence 
of the Lagrangian with respect to the gauge potential and its derivatives. 
To solve this hierarchy is the main aim of Utiyama's general program. The 
solution of the above system gives the covariant objects of the theory as well 
as the functional dependence of the Lagrangian on these objects. 



3. 1 Hierarchical Equations Solution 



3.1.1 Solution of the Equation (f77|) 



Once equation (fT7j) involves only the dependence of the Lagrangian density 
on the second derivatives of the gauge potential, one may propose that this 
dependence appears through an object R a QUti which must have a cyclic per- 
mutation symmetry: 

^ + + R a m v = 0- (18) 
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The most general linear object constructed from d 2 A with this property is 
given by: 

R a aufl = d v d a A^-d a d,A a v , (19) 
where we chose A = 1 without loss of generality. 



3.1.2 Solution of the Equation (G] 
Let us consider now the equation (fT6|) written in terms of (A, dA, R): 



dL 



(i) 



d (d u A c 

d(d a A\ 



+ 



+ 



dL 



(i) 
o 



04 



(i) 



d[R a 
dL { 



(1) 



d [R a , 



avfi 



d (R% u ) 
d (R% a ) . 



f c \A\ = 0. 



This equation shows a symmetry in v <-> a, which must be present in its solu- 
tion. Also, the solution will be construct from R and then it must respect the 
cyclic symmetry as dieted by (ITT)) . Therefore, the solution for this functional 
partial differential equation is such that the functional dependence of Lq with 
respect to R shall be through the object 



Q d f3 P o = R d a P <j ~ feb A a d p A°0 — A b d a A 



(20) 



With this new quantity we pass from Lq (A, OA, R) to L Z> (A, OA, Q). 

A second independent solution to (TTBT) can be found if one rewrites it in terms 
of 4 2) (A,dA,Q), i.e., 



'(2) 



dL& 



dL 



(2) 



d(d v A° a ) d{d a A% 



= 0. 



In this case, the dependence on the second derivative is eliminated and the 
dependence on the first derivative can figure only via the combination 



Al 



d u A d a A 



Thus, when we construct Q we are actually selecting a sector of the gauge 
potential Lagrangian. Here we have introduced the antisymmetric operation 

( t }\ 

O^u,] = 0^ — O vil for a general 0^ u . So, one goes from L y Q> (A, dA, Q) to 



r(3) _ r(3) / m ac ^ 



(3pcr 
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3.1.3 Solution of the Equation 173]) 



Our next step is to rewrite the equation ([15]) in terms of the new functional 
dependence of the Lagrangian, L = (^A, A[],Q^: 



dL 



(3) 



+ 



dL 



(3) 



°— (S" p S a a - 6^ f a d e A e a + 



(3) 



+ 



We*, 



( S a A \pf3] + ^ A \pa] + S p A ya}) fac + A % {^a^p ~ ^ pK) f/cfa, 



= 0. 



The solution of this equation utters that the functional form of the Lagrangian 
density depends on three objects: 

F d pa = Af pa] + f a d e A e p A a a , (21) 



G% pc = Q% a - {{^A c m + S^A c [pa] + 5iA\ M ) f d c A\+ 

- (^p ~ A \^ f b d J g C e A \ Ae p} (22) 

and A a „ itself. 

Keeping in mind Lq^ = {F^ ■> ^%pa^ we reexpress (TT51) : 

dL^ (F, G,A) _ 
dA a v 

Due to this, the Lagrangian kernel cannot be explicitly dependent on the gauge 
potential A, which is the real reason why gauge fields are massless. Therefore, 
the presence of massive terms is only possible if the gauge symmetry is broken. 



3.1.4 Solution of the Equation ( f77| ) - Condition on the Lagrangian 

After all that, we must put the equation (TTJJ in terms of the new objects F 
and G, 

d Fr fsb [dA\ A » + d ( dvA tf" A <) + (23) 

9Lq_ f h\ dGd (3pv Ab , dG %fxr a Ab , dGd f3pa a ~ A b 

+ dG^ L b [-dA^ A " + ^ < + d (d idv A^ A C 
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or, applying the Leibniz rule and the Jacobi identity, 

fll-(4) or (4) 

This is a condition on the Lagrangian of the gauge potential which must be 
satisfied for the local symmetry to hold. 



3.2 New Expression for G and the Bianchi Identity 



In order to reduce the number of independent objects in the expressions above 
we will write G in terms of F and A. To do so, let us consider initially that, 
from (USD and (J2H), we get 

R h p pa = dpF% + f e \dp (A e p A b a ) . 



Using these informations to evaluate G it results 



GV = d p F% + f c \ [A%d p A h a + A\d a A c p\ + 

+ fJb (^ 6 fl^'[po-] + ^V^W]) ~ fgbfcf-^p (A b a A c p — A h p A c a 
that, with (1211) and the Jacobi identity, conduct us to 



G a 



Ppa 



9pF% + f c a e A e 3 F c oa = D%F 



c e P pa 



cP pai 



where Dp is a kind of Fock-Ivanenko derivative [12 



c/3 



(25) 



(26) 



This kind of structure in terms of a Fock-Ivanenko-like derivative is exhibited 
also by F when written as 

F% = D d ap A\ - D d aa A a p = D d a[p A a a] , (27) 

in which 

D d ap = 5 d a d p -u d ap - u d ap = \f d a A%, (28) 

with a factor 1/2 in the spin connection uj d p , different from eq. ([2"B"j) . This 
difference refers to the fact that A e is a connection and not a vector like F d . 

With this new expression for G we immediately verify that it is antissymet- 
ric in its last two spacetime indices. This is consistent with the antissymetry 
property required by equation (fl6j) . On the other hand, equation (|T7|) stab- 
lishes that a cyclic property must be present in those objects containing the 
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second derivative of the gauge field. In fact, starting from (|25l) and with the 
expression of F d a and the Jacobi identity we find 

GV + G V/3 + G %p = + D % F lp + D a ca Fl p = 0, (29) 

the so called Bianchi identity. 

Thereafter, the Bianchi identity ff29l) is a natural consequence of the local sym- 
metry in the second order formalism as dictated by the hierarchical equations 
(|14p to (|17p . and not just an a priori equality constructed with F. 



This identity and the manifest skew-symmetry of F a in its spacetime indices, 
allow us to rewrite G as 

G%kt = D^Fpyj, (30) 
in the same suggestive form as (T27j) for F. This kind of structure will be 
explored below in the section 16. 1[ 



3.3 Transformation Laws 

Next, we construct the transformation law for G. Remembering [1], we have 

&A\ = fc a b A b ,e%x) + — 

and 

SF% = e c (x)f« b F%. (31) 

Therefore 

*GV = e d f d a g d p F^ + e d f c a J d c g A e p F% + f d a g F%d p e d + f \SA%F% 

and, applying the transformation law for the gauge potential and the Jacobi 
identity, 

5G% a = e d (x)f d a c G% a . (32) 

Then just like F, G is also covariant under the action of the local Lie group0 
i.e., it transforms like a vector under the action of the Lie group. 



(4) 



1 We could have arrived at the same conclusion with a glance at ([23]). 9F ° d f d hFp\ + 
-qq£ — fc h Ppt ^ 0' wn ich is nothing but the invariance condition of the Lagrangian 

fjpa 

8Lf = 0, %^8F% + ^~JG% a = 0, and so, 5F% = e c f d h F^; bG% pu = 

^ C 3 c}fi h pp(,i where we have used the fact that the Lagrangian is a kernel functional 
of F and G alone. 
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4 Second Order Current 



From the minimal coupling prescription the total Lagrangian is 

L T (Q A , V m Q a , A, dA, d 2 A) = L (Q A , V,Q a ) + L (F, G) . 

Observe that the total Lagrangian is a function of the second derivatives of 
A, by means of G in L . 

From the variational calculus, 

5L 8Lt „ 3Lt 
W = dQ X ~ u d(d u Q A Y 

5L T _ dL T I dL T \ ( dL T \ 

The vanishing of the variation of the total Lagrangian results in an equation 
composed of two terms: a volumetric and a surface one, which must be null 
independently, in a development similar to Noether theorem. 



4-1 The Variation of the Total Lagrangian 



The variation of L T = L T (Q A , dQ A , A, dA, d 2 A) is: 



SLr = 



5Aa uA,e ^y 5A% 



e a + 



5L 



r 



6Q< 



5Q J 



+ d u M\ 



where 



M v = ^> + ^^6Q A 



5A\ d{V u Q A ) 



1 f d Ab I ~~U , 

Jab p \ dG d ,, . , dG d 



dL 



(4) 



dF 



[Vfl] 



dL 



(4) 



+ 



dL 



(4) 



dL 



(4) 



+ 



dL 



p[vp] 
(4) 



SAI+ 



d P (SA^) , 



or 



M u = N c »e c + Ord,e a + -P a ^d p d,e a , 
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with 



5Ln 



+ 



+ 



5A% 1 d(V u Q A , 

ara ~ Jab p 



-T, 



(c) B<2 + 



+ 



"5* 



(1) 



+ 



v[pp] 



dG a „ 



" 2 



(4) 



+ 



«9L 



(i) 
o 



v\pp\ 



OG 



(33) 



pW 



b . 

A 4 ' 



(4) 



QUp 



dF lM 



2*' 



<9G a 



df 
dG %iA 



dG %iA , 



+ 



dG a 



p\vp\ 



1 

+ 2 



<9L 



(i) 



<9G 



i/[/xp] 



dG 



H[vp] 



(34) 



fab^p ) 



and 



(4) 



+ 



^4 

dG a , 



(4) 



(35) 



^.S Hierarchical Equations for the Current 



The vanishing of the volumetric term in 5L? = gives the equations of motion 

6Q I = 0, SA^ = ' 

By the same way from the vanishing of the surface term, 

d v M v = 0, (36) 

we have the following set of hierarchical equations: 

d u N c u = (37) 

N e " + d^O/" = (38) 

O a ( ^ } + d p P a u/ "* = (39) 

p a pUfl + p a w + p / pu = (40) 



These equations constitute a hierarchical set of equations which governs the 
conservation law associated to the local symmetry, in a similar way that eqs. 
(j!4j) to (j!7p determined the functional form of the Lagrangian L . 
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4-3 Solution of the Hierarchical Equations for the Current 



Equations (flOl) and (I39I) are automatically satisfied by virtue of the symmetries 
in F° v and G\ pu . 

Condition will be used to defflne the current, while the equation (13T1) sets 
the conservation law. 



4-3.1 Qu&siconserved Current 
Let us define, as done by Utiyama 
8Lt dL T 



J' 



dA c 



+ D a , 



d(V u Q^ BQB 



dL 



(i) 



flr(4) 
(4) 



dGi 



_ fd rpb _|_ 

J cb pu' 



upp 



dp dG 



(41) 



This definition is inspired by the most direct experimental sense of current, 
as a measure of the response of the system under a variation of the field. 



The condition (1371) now gives 



d v J c l 



d ^- l _ d d 
5A C „ 2 v p 



dG d 



dG d 



p[vp] 



f d A h 

lcb A p 



(42) 



On the mass shell it follows the quasi conservation of the current J£. The term 
"quasi" is understood in the sense that if one takes the integral of f|42l and 
chooses boundary conditions such that 



,-><i 



da v d p 



dL 



(i) 



+ 



9L^ 



dG 



rd Ab 
hzb A p 



the conservation of is globally recovered. For instance, this kind of bound- 
ary condition occurs when A e p and its first derivative are null on the boundary 

da 



It is interesting to notice that in the case of an abelian group this conservation 
is also achieved. The same occurs in the first order approach, and that is 
because Utiyama's definition coincides with Noether's current. This does not 
happen in the second order theory, as one can immediately see from the non- 
conservation of J£. In order to establish a conserved current in the second 
order approach, we must build another one based on Utiyama's proposal, 
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which is later compared with the standard Noether's current. This is done in 
the following section. 



4-3.2 Conserved Current 

An alternative definition for the current is 



V 



dL 



T 



dA c 



-d 



dLn 



d{V v Q A 



(c) bQ" — 



9Lt _ t a r\U 



dG d 



f <* TP" _|_ 
J cb pp' 



vpp 



I D d 

1 U a p dG d 



(4) 



p[up] 



fcb^p ^P 



ar(4) 0.(4) 
°' L , r d Ab °- L 



(43) 



Therefore 



d„J r v = d v 



5L 7 
~5A C 







under the equations of motion, showing that the current (143]) is strictly con- 
served. 



4-4 Concerning the Utiyama and Noether's Currents 



In the current ( HIT) the first two terms are the current obtained by Utiyama in 
the first order formalism pQ, the third and fourth terms are of second order, 
and the last one is the gimszconservation term. 



On the other hand the current (j4*3|) is composed by first two terms of the first 
order current, plus second order terms involving G, and the last is a topological 
term [13] . in the sense that it is conserved independently of the equations of 
motion. 



It is worth to remember that these topological terms cannot have their origins 
explained by the dynamics, but the general local invariance scheme brings 
them with it. 



We can evaluate the transformation laws for fill and fl43l: 



5J„ 



e fec^a 



34 
dF iM 



d a h 



+ 



pWp] 



dG d . 



p[up] 



+ OL ° f a d d e e 
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dL 



(4) 



dL, 



(4) 



f a d e 



dL, 



(4) 



- f a d d e e 



f)(~<a Jec^^P^ 



As expected from J» = and J" = - ^ a( gj eJ 
not covariant. However, one can define a covariant current: 

dL (Q, VQ) <9L T 



both currents are 



-T C A B Q B , 



dA^ d (V P Q A ) 

which is not strictly conserved, but only covariantly conserved: 

D a j M = 0. 



(44) 



(45) 



We conclude that covariance and conservation never can be obtained simulta- 
neously; in order to maintain one we must sacrifice the other. 

To compare the proposed currents J" and J" with Noether's one, we rewrite 
them as 



7 " 



dL 



T 



t A n B - 



dL n 



+ 



d 



dL 



d (d u d p A^) 
dL 



f c a b d p A b ^+ 



P d(d p d u A^) d(d v A%) 



dL 



and 



7 v 



dL 



T 



T, A " B 



+ 



d, 



did„(r)- 

dL 



Q - 



dL, 



d (d u d p A^) 
dL, 



fc a b d p A b 



dL ri 



'd(d p d v A^) d(d u A^) 

which are directly comparable with the Noether's second order current, 

dLf ^ a 



(Jn) c 



t A n B - 
d(d u Q A ) 2 ^ BQ 



+ 



8r. 



dL 



d (d u d p A^) 
dL 



f c \d P A\+ 



d (d p d v A^ d (d v A^) 



Jc b A is 



Then, one easily sees that the expressions of J and Jn differ between them 
only by the presence of topological terms, while in J there is a quasiconser- 
vation term, d u 7 9L ° — r-, which does not contribute in the first order case. 

Therefore, in the first order approach, Utiyama and Noether's currents coin- 
cide, but in the second order case it is necessary to perform a generalization in 
order to accomplish the conservation, giving rise to the presence of topological 
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currents. These topological factors are naturally introduced by the presence 
of the additional term fyt g/g 1 ^ ) i n J e v > which can be interpreted as a flux 
of gauge potential, extending the usual definition of current as a variation of 
"energy" L T with respect to the external field A c u , J^r- 

Another way to see how to implement this generalization is to remember that 
the conservation law of the current is intimately related to the last of the 
hierarchical equations. For instance, in the first order approach, 



where we have used the equations of motion, and from where follows the 
conservation as a direct consequence of the equation 

dL dLp _ Q 



Analogously, in the second order case we extract from the equation of motion, 

= „ dLr_ _ a dL dL 
Jc dA% »d(d„A%) ° p0fl d(d p d,A%y 

which again is conserved by virtue of the equation (fl7|) . 

On the other hand, Noether's procedure for the calculation of the current 
imply exactly the first of our hierarchical equations, (1371) . giving the Noether's 
current in terms of F and G as just the N" — fjf- quantity. We observe 
again that this, up to topological terms and a global signal, is the conserved 
current J". Therefore, we see that the Utiyama's systematic method leads 
to the expected result and gives the bonus of finding topological currents. 
We conclude that far from being arbitrary, this topological current is induced 
by the structure of the hierarchical equations, arising from the local gauge 
invariance principle and has direct observable consequences for the charges. 
This will be illustrated in the next section. 



5 Applications 



5.1 U (1) Group: Podolsky's Generalized Electrodynamics 



Usually in the literature [9] Podolsky's electrodynamics is supposed to be 
the simplest generalization of Maxwell theory whose Lagrangian, contain- 
ing second-order derivatives of the electromagnetic potentials, is gauge and 
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Lorentz invariant and still leads to linear local field equations. In fact, here we 
will show that the second order gauge theory is able to prove that Podolsky's 
electrodynamics is the unique theory which has such properties. 

Podolsky's second order theory for the electromagnetism consider the La- 
grangian density [8] 

L = -\f, v F» v + jdTF^dpFi*. (46) 



We will check that the Podolsky's theory fulfils the condition for a second order 
gauge theory. The equations to be satisfied are ([14"]) to ffTTj) . The first equation 
is automatically verified for the U (1) group by virtue of the nullification of 
the structure constants. Therefrom, for this group any scalar constructed with 
the tensors F and G is a possible Lagrangian, in principle. 

The next two equations are identical to the first order ones 

dL 



3A h 

dL dL 



0. 

= 0, 



d{d v A a ) d(d a A, 
and consenquently hold in the case of Podolsky Lagrangian ( 1461) . 
It remains to be verified the last of the hierarchical equations. From ( 146 

, , = aY" (cfcFA e - (Fd £ A°) - r <- a, 

d {d p d a A T ) 

and by cyclic permutation, 

dL dL dL 



d(d p d a A T ) d(d a d T A p ) d(d T d p A a ) 

and we see that the Podolsky's theory is in fact a second order gauge theory 
a la Utiyama. 

Alternatively, one can start straight ahead from the second order gauge theory 
for the U (1) group and write the G tensor (|25|) . 

G/3 pa = dpF pai 

from which we can construct a vector, e.g., 

G% = d^Fpp = G P 
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and obtain a second order Lagrangian equivalent to the additional term pro- 
posed by Podolsky: 

Notwithstanding, there are other possible Lagrangians densities. For instance, 
considei[!] 

L G = b 2 G^G Ppa = b 2 ^F^d p F pa , 

which satisfies all the hierarchical equations by the same arguments given 
before. 

Using the Bianchi identity (1291) one concludes 

G^G Ppa = 2d^F^d a F pP . 
Perceive that this is the Podolsky Lagrangian apart from a surface term: 
d^F^d a F pP = d a (d^F^F^) - (dpF^Fpp) + d a F^d^F pP . 

Then L G is equivalent to Podolsky's taking 2b 2 = 

We can explore another order of contraction to the indices of G: 

L G = c 2 G^G pPa = c 2 d^F^d p F aP . 

By the cyclicity symmetry fl29l) 

L G = c 2 d p F^d p F^. 

Again, this is Podolsky Lagrangian, except by a surface term, with c 2 = ^L. 

Therefore we can infer that all the quadratic Lagrangian in G are reducible 
to the Podolsky's form since there are only two possible cases for indices con- 
tractions with respect to the derivative of F: or the indices in the derivatives 
are contracted between them or they are contracted with one of the FF. In 
the first case one uses the Bianchi identity, which reduce all contractions to 
the second case. But the second case is just Podolsky's, or can be put in this 
form by an adequate surface term. 

For these reasons, we see that the most general Lagrangian for the U (1) group 



2 Of course the complete Lagrangian for the gauge field is the Maxwell one plus the 
term under discussion. 
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in a theory a la Utiyama is 

which is the Podolsky Electrodynamics Lagrangian. 

This proves that Podolsky Lagrangian is the unique linear second order gen- 
eralization from Maxwell theory compatible with the gauge principle. 



5.2 U (1) Currents and the Mass of the Photon 
With flUJ) for the case U (1), 

jv _ -v _ dL T T A B 

and this means that the second order formalism has no effect on the current 
if the symmetry is U (1). 



Now if we employ (j43 



8T 

r = f - <9 M -— ^- = f - d^, 

which differs from J u by a topological term. 

The equation of motion, in vacuum, for Podolsky's electrodynamics is 

(l + a 2 D) <9 A i^ A (x) = 0, 
which, with the Lorentz gauge condition, d^A 11 = 0, reduces to 

(l + a 2 D) DA" (x) = 0. (47) 
Extracting the Fourier transform, we find two possible dispersion relations: 

p 2 = 0, 

p 2 -- = 0. 

a z 



The first one corresponds to a massless mode, E 2 — p 2 = 0, while the second 
is a massive solution: 

1 I 

~7i 



E 2 = p 2 + ^ , m = -. (4J 
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In the case of a static isolated point charged particle, equation (T4T|) gives: 

-(l-a 2 V 2 ) VV (r) =47rp(r), 
p (r) = — e5 (r) . 

We have selected a Lorentz gauge such that = (<p, 0), and used the signa- 
ture f] = diag (1,-1,-1,-1). Via Fourier transform, we find 

e / _r 

V ( r ) = — (1-e « 

The finite reach of the massive mode suggests the presence of a shielded region 
around the point particle. 

The conserved current is given by 

j°(r) = -e(5(r) + -±-e-S) , J k = 0. 
V err / 

Integrating this density over the whole space, we obtain the conserved charge, 

q = -e(l- Ait) 

which is just a simple renormalization of the naked charge. However, follow- 
ing Moniz and Sharp [TJ], the quantum theory of the nonrelativistic electron 
provides a natural cutoff at short distances whose magnitude is of the order of 
the Compton length e. Employing this idea of a shielded region, the effective 
conserved charge obtained is 



leff = -e 1 - 4tt 



a + e ( e 

exp — 

a \ a 



(49) 



This effective charge can be applied to estimate the mass of the photon, ac- 
cording to (HHI) . Equation (j4"9"l) can be rewritten as 



jrd - leff -e _ a + e ( e 



a™ = ^ = Air— exp 



This quantity can be interpreted as the relative uncertainty of the electron 
charge. Solving this transcendental equation for Podolsky's parameter a, we 
find a measure of mass for the massive photon, 

h 

171 = . 

ac 

With experimental data from [15] and a simple propagation of errors neglecting 
the uncertainty in (J r e el , we obtain 

a = 1.105 868 617 (14) x 10~ n cm , 
m = 1.985 370 21 (17) x 10" 13 eV . 
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This estimative cannot be ruled out in accordance with experimental data 
reported in [T5] and [16] considering null tests of Coulomb law. These refer- 
ences give m < (9.0 ±8.1) x 10~ 10 eV. In the context of an effective classical 
theory, this mass must be interpreted as an energy scale characterizing the 
regime where the Podolsky's electrodynamics becomes relevant, in a similar 
way as the phonon mass characterizes the effects of particle excitations in a 
cristal [17] . It is meaningless to think of a rest mass for a phonon, since it is 
an excitation of the vibrations in a cristal. The same interpretation can be 
applied for the photon in Podolsky's theory. 

We emphasize that these results were derived using a completely classical 
approach and not following a construction of effective classical theories from 
their quantum versions. This is the reason why the renormalized charge is 
found and no estimative for the photon mass can be made in the absence 
of a shielded region. Even taking a shielded region, if Noether's current was 
considered, the effective charge would be just e and, again, no mass for the 
photon could be found, as well no renormalization for the charge would be 
obtained. 



5.3 SU (N) Group: The Lagrangian AAB 



As a non-abelian example we apply the second order gauge theory to the 
description of Alekseev-Arbuzov-Baikov's effective Lagrangian [18] proposed 
to eliminate infrared divergences in SU(N) theories. In our notation, it is 

i i 

t /~ia r~i \uv i e a rpbXp 7-1 r 771c 

where M is the mass scale of the infrared gluon. This scalar density satisfies 
the condition (j24p . a fact verified after using the total antisymmetry of the 
structure constants. 



The conserved current for this theory can be obtained from (1431) : 

{^)>wM*'-'~ 1 



7-j d s-i pup 1 f d rpepr rp 1 
a p^d ' Jae dr 



2 pn a I 



j ^_f d f) f A b (1 pufJ/ -\- F btMT F u 

+ Jc b°P p^d + r r dr 



which differs from the current in [TT] only by topological terms. 

The equivalence of 1-loop infrared regime and the second order gauge theory 
shows that an analysis of high order gauge theories could give more information 
about higher loop expansions. 
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6 General Remarks 



In this work we saw how Utiyama's approach to gauge theory can be extended 
to include second order derivatives of the gauge potential. The development 
showed the need for the introduction of a new object G c p pa , as well as the 
already familiar field strength F% a . Both are covariant and have a Bianchi 
identity associated with it, which was derived from hierarchical equations for 
the auxiliary field Lagrangian. 

The nullification of the variation of the total Lagrangian conducted us to 
conserved currents that have two possible definitions. The first one was defined 
as in the first order treatment, but it is only quasiconserved. Conversely, we 
could define a new current, which differs structurally from the one proposed by 
Utiyama, but is characterized by its conservation. Another advantage of the 
general Utiyama's procedure is the obtainment of topological currents from a 
local invariance. It is a new tool to study topological aspects of gauge theories 
from their local symmetry, and must be better investigated in the future. 
These topological currents are not arbitrary, but rather implied by the very 
structure of the field and hierarchical equations, as will be better discussed in 
the section 16. 1[ 

We also saw how the second order development can be implemented in the 
case of the U (1) group leading to the well-known Podolsky Generalized Elec- 
trodynamics, whose Lagrangian was demonstrated to be the unique possible 
extension up to surface terms. As an example for the more complicated case 
of a non-abelian theory, we applied the second order gauge theory to the 
Alekssev-Arbuzov-Baikov effective Lagrangian valid in the infrared regime, 
demonstrating that the conserved current (1431) and the one obtained in [TT] 
differ by surface terms. In this way, we have shown that 1-loop effective terms 
in the action can be found from first-principle calculations imposing the gauge 
symmetry to second order Lagrangians. In addition, an extension of Utiyama's 
method to higher orders derivatives could give important information about 
higher loop expansions. 

The conserved current laws were used to estimate the energy scale of the pho- 
ton massive mode, in the electrostatic regime, and the result is in accordance 
with the present observed limits. This massive mode was engendered by the 
topological terms in the conserved current (|43|) in accordance with the known 
fact that topological terms are responsible for dynamic mass generation mech- 
anisms [19]. Besides, the conserved current J gives the renormalized charge, 
which is quite natural in view of the equivalence among second order gauge 
theory and 1-loop expansions. On the other hand, the Noether's current only 
gives bare charges. Notice that in the U (1) case both J and J are conserved. 
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Among the perspectives opened by this work we can cite another natural 
extension to Utiyama's work allowing the Lagrangian to depend also on the 
second derivative of the matter field d 2 Q, which is presently under study. 
Other possibilities are commented in the next sections. 



6.1 Generalization for Higher Order Theories 



In order to get a generalization to higher orders, we note that some structures 
apparently repeat themselves in the first and second order formalisms. For 
example, the introduction of the gauge potential as a compensator to achieve 
the local invariance and the hypothesis L = L (A, dA) leads to conclude that 
the derivative of A must respect the symmetry 

dL dLn , . 

+ —, — = 0. (50) 



d(d,Al) d [d vA 

In the second order extension one finds that the second derivative accomplishes 
dL ° + + ; - 0. (51) 



Therefore, it seems natural to suppose for a nth-order theory the higher order 
derivative appearing only in an object 

pa = Q{ai,...,a„}Q a Aa 

lh a 1 ...a n — u Oil ■■■ u «n-i n a„ 

P{ai,,.,,a n } 

carrying a cyclic permutation symmetry, 

H Ra ai ...a n = 0, 

P cy cl{ai,---,a„} 

where P a n } is to denote the permutation of indices and P cyc i «n} 

the cyclic one. This identity might be understood as a generalization of the 
Bianchi identity. 

Other recurrent structures are as follows. Introduced the auxiliary field A, the 
tensor F can be written as (I2T|) 

In the second order case a new object comes into place (I3U1) 

Go r-ia Z? c 
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In time, we may conjecture for a third order theory the emerging of a new 
object whose functional form is 

n aj3pa — u c]pS* fi\pai 

respecting, as said before, a Bianchi identity 
Remembering that for a general Lie group vector, 

7-) a r* c [n n l a T?d fa j^e rpd 

U c[p^ 0\pa - L /*> U fi\ d P pa ~ Jde r pfi? pa, 

one finds 

= ^%p C p}fm + ^ a c [/3^ C p]aiJ, + ^ >a c[p ( ^ C <7]pf3 + -^ ) c[ t r ( ^ C / i ]/3p = 

X a pe , y a pe j a pe r a jpd ne pi 

— J d e r pf3 r pa ~1~ J d e r /3p r ap J d e r pa r p/3 J d e r ap r f3p ~ u • 

and the structure is in fact consistent. It is important to note that this is not 
a rigorous proof of the third order structure, but a significant indication of its 
form. 

We have seen that the local gauge invariance condition on the first order 
Lagrangian is [1] 

;)/ 

fdph — n 
Jch P pa — U 



dL 



f)pd J cH pa 

and we have proved for the second order theory the generalization 

9L rd F h , 9L Q d h _ 
dF% Jch pa dG% p J c h Pp ° ~ 

Therefore, it appears natural to expect the restriction 

dL o f d F h , 9L d h OLp f d H h _ n 

flpd Jch r pa ' p\(~id Jc h^" f3pa ' a rrd Jc h n a/3pa — u 
OP pa OLr f3pa OU a/3pa 

on the third order Lagrangian of the auxiliary field. 

The structures identified until now indicate that for still higher order theories 
one can hope for objects given in terms of Fock-Ivanenko derivatives of the 
quantities of the former order. Another feature is that the higher order objects 
also are supposed to satisfy a Bianchi-like identity. 

Regarding conservation laws, Utiyama defined the current 



J? 
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which is a conserved quantity for the first order theory, Lq = Lq (A, dA), under 
the equation of motion 

8L _ dL T _ dL T _ 
6AI ~ dA^ v d(d v Al) ~ 

We have seen that this definition of current is not appropriate for the second 
order approach, L = L (A,dA, d 2 A), since it is not a conserved object. To 
solve this problem we defined 

a ~9K V d(d„A$ 

as the expression of the conserved current valid on mass shell, which in this 
case reads 

5L dL T dL T dL T 

-n— - TTi °v — 7 r + OpOu — 7 r = 0. 

8AI M« d(d u A^) P d(d p d u Al) 

Notice that in each case we have just isolated the last term in the Euler- 
Lagrange equation and defined the remaining terms as the conserved current. 
The conservation is, in fact, assured by the symmetry imposed by the last of 
the hierarchical equations in each case, namely (150]) and (15T|) respectively. 

This reasoning may be applied to extended theories: in a n-order theory Lq = 
Lq (A, dA, d n A) the current may be defined as the first n terms of the Euler- 
Lagrange equation. The currents constructed following this recipe would differ 
from the Noether's one by topological terms, a fact that reflects upon the 
values of the charges. Which values is the correct one stands to be selected by 
experience. 

6.2 Geometrical Aspects 

Geometrical aspects are not evident in the Utiyama's approach since it is 
essentially an algebraic implementation of the local symmetry principle. A 
geometrical interpretation of G could help one to understand the recurrent 
emergence of the Fock-Ivanenko derivative and the general structure of higher 
orders Lagrangians outlined above. 

Moreover, the comprehension of the geometry seems to be fundamental if one 
wants to describe higher order gravitation as a direct application of Utiyama's 
second order procedure, a work now under construction by the authors. This 
could also illuminate some features of the nonrenormalizability problem in 
gravitation. 
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6.3 Application to QCD 

Besides the analysis of the AAB Lagrangian, it would be interesting to per- 
form an 1-loop calculation and check if the cancelling present in Podolsky 
Electrodynamics [8] is also manifested in QCD. 

However, one must consider which quadratic form in G will be explored. Due to 
the nonlinearity of the theory, it is possible that all combinations be necessary 
in order to obtain a renormalizable theory. Another probable feature is the 
impossibility of closing the number of derivatives (counterterms), as in the case 
of gravitational field. Some tentative works in this direction were undertaken 
in [TO- 

An additional point is to investigate if higher order terms gives more informa- 
tion about the confinement phenomenon using, for example, Wilson criterion, 
as noted in [IT] . Once more, to recognize the geometrical connection which 
defines the holonomies would be a good guide. 



6.4 Constraint Analysis 

The local gauge symmetry imply the existence of constraints in any order of 
derivatives which require a special attention in order to construct the Hamil- 
tonian description. 

The analysis of the constraints for a second order Lagrangian can be imple- 
mented using several approaches, such as that in [20] for Podolsky Electro- 
dynamics. Nevertheless, the study of an specific Lagrangian for non-abelian 
groups is interesting in view of the applications for QCD and gravitation. 
This study is presently under implementation by means of Hamilton- Jacobi 
technique [21]. 
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